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Abstract 

A  two-phase  flow  model  for  a  proton  exchange  membrane  (PEM)  fuel  cell  is  developed.  The  model  is  based  on  the  mixture  flow  model 
and  the  unified  approach  is  used.  Instead  of  using  a  separate  model  for  the  catalyst  layer,  the  catalyst  layers  are  now  included  in  the  respective 
unified  domains  for  the  cathode  and  anode,  thus  continuity  boundary  conditions  at  the  interface  between  the  catalyst  layer  (CL)  and  the  gas 
diffuser  layer  (GDL)  are  no  longer  needed.  The  model  couples  the  flows,  species,  electrical  potential,  and  current  density  distributions  in 
the  cathode  and  anode  fluid  channels,  gas  diffusers,  catalyst  layers  and  membrane  respectively.  Furthermore,  the  two-phase  flow  model  is 
also  used  in  the  anode  side,  and  the  momentum  transfer  between  the  liquid  and  gas  phases  due  to  phase  change  is  taken  into  consideration. 
Experiments  have  been  conducted  to  study  the  performances  of  a  PEM  fuel  cell  and  the  results  are  used  to  improve  and  validate  our  model. 
The  modeling  results  of  polarization  curves  compared  well  with  the  experimental  data.  The  model  is  used  to  study  the  influences  of  fuel  cell 
operating  temperature,  operating  pressure  and  humidification  temperature  on  the  oxygen,  vapor  and  liquid  water  transports,  as  well  as  fuel 
cell  performances. 

©  2004  Elsevier  B.Y.  All  rights  reserved. 
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1.  Introduction 

It  is  well  known  that  water  management  is  a  critical  issue 
in  proton  exchange  membrane  (PEM)  fuel  cell  design  and  op¬ 
erations.  On  one  hand,  the  polymer  membrane  in  PEM  fuel 
cells  needs  to  be  hydrated  to  maintain  high  proton  conduc¬ 
tivity.  On  the  other  hand,  excessive  liquid  water  may  occupy 
a  significant  fraction  of  the  pore  volume  in  the  catalyst  layer 
and  gas  diffusion  layer  (GDL),  thus  may  significantly  reduce 
reactant  transfer  rate  to  the  reaction  sites  and  cause  flooding. 
Therefore,  proper  water  balance  and  management  is  essen¬ 
tial  to  maintaining  high  performance  of  a  PEM  fuel  cell.  A 
good  understanding  and  proper  modeling  of  the  various  trans¬ 
port  processes,  especially  water  transport  is  essential  for  the 
improvement  and  optimization  of  a  fuel  cell  and  fuel  cell  sys- 
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terns.  At  present  many  of  the  transport  phenomena  inside  a 
fuel  cell  cannot  be  directly  observed  or  measured;  thus  mak¬ 
ing  mathematical  modeling  a  critical  tool  to  understanding 
these  transport  phenomena. 

Some  of  the  pioneering  work  in  PEM  fuel  cell  modeling 
includes  those  by  Bernardi  and  Verbrugge  [1,2],  Hill  and  Ver- 
brugge  [3],  Springer  et  al.  [4,5].  These  one-dimensional  mod¬ 
els  laid  the  foundation  for  PEM  fuel  cell  modeling.  Pseudo 
two-dimensional  model  includes  those  by  Fuller  and  New¬ 
man  [6]  and  Nguyen  and  White  [7],  which  included  the  ef¬ 
fects  of  reactant  consumption  along  the  flow  channels.  Gu- 
rau  et  al.  [8]  developed  a  two-dimensional  PEM  fuel  cell 
model  that  coupled  fluid  flow,  mass  transfer  and  the  electro¬ 
kinetics,  thus  introducing  the  computational  fluid  dynamics 
(CFD)  into  fuel  cell  modeling.  The  approach  adopted  was  a 
unified  approach  that  eliminated  the  need  of  prescribing  arti¬ 
ficial  boundary  conditions  at  the  various  interfaces  between 
the  different  components  of  a  fuel  cell.  Zhou  and  Liu  [9] 
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Nomenclature 

a  catalyst  surface  area  per  unit  volume 

(cm2  cm-3) 

C  species  mass  fraction 

D  diffusion  coefficient  (cm2  s-1) 

E  potential  (V) 

Em  Equivalent  weight  of  ionomer  (g  (equiv.)  - 1 ) 

/e  surface  area  of  liquid  water  per  unit  volume 

F  Faraday  constant  (96487  C  mol-2) 
io  exchange  current  density  (A  cm-2) 

I  current  density  (A  cm-2) 

J(s)  capillary  pressure  function 

kr  relative  permeability 

K  absolute  permeability  (cm2) 

L  characteristic  size  (cm) 

M  molecular  weight  (kg  mol - 1 ) 

N  mole  flux  (kg  cm-2  s- 1 ) 

p  pressure  (Pa) 

R  universal  gas  constant  (8.314  J mol-1  K-1) 
s  phase  saturation 

Sh  Sherwood  number 

T  temperature 

u  velocity  vector  (cm  s  - 1 ) 

u  the  v  direction  velocity  component  (cm  s  - 1 ) 

v  the  y  direction  velocity  component  (cm  s  - 1 ) 

v  mole  fraction 


Greek  symbols 

a  net  water  transport  coefficient 

aa  anodic  transfer  coefficient 

oic  cathodic  transfer  coefficient 

y  multiphase  correction  factor 

8  width  (mm) 

s  porosity 

t)  overpotential  (V) 

X  water  content  in  membrane  or  individual  mo¬ 

bility 

li  viscosity  (N  s  cm-2) 

v  kinematic  viscosity  (cm2  s-1) 

p  density  (kg  m- 3 ) 

cr  interfacial  tension  (N  cm-1) 

Subscripts  and  superscripts 
a  anode  or  air 

act  activity 

c  cathode  or  capillary 

d  electro-osmotic  drag 

diff  diffusion 

eff  effective 

g  gas  phase 

hyd  hydraulic 

in  inlet 

k  phase  k 


1  liquid  phase 

m  membrane  or  mixture 

r  relative 

s  solid 

sat  saturation 

v  vapor 

w  water 

a  species 


and  Um  and  Wang  [10]  developed  three-dimensional  models 
using  a  similar  approach  by  Gurau  et  al.  [8]. 

Wang  et  al.  [11]  and  You  and  Liu  [12]  developed  PEM 
fuel  cell  models  that  applied  two-phase  flow  model  for  the 
cathode  side  of  the  PEM  fuel  cell.  The  model  by  Wang  et 
al.  [11]  treated  the  cathode  catalyst  layer  as  an  surface  and 
did  not  included  the  membrane  and  the  anode  side;  while 
the  model  by  You  and  Liu  [12]  treated  the  cathode  catalyst 
layer  separately  by  a  pseudo-homogenous  model  [13]  and 
only  single-phase  flow  model  was  used  at  the  anode  side. 

In  this  study,  we  developed  a  two-phase  flow  model  using 
a  similar  treatment  as  in  You  and  Liu  [12]  and  Wang  et  al. 
[1 1],  but  our  model  treated  the  catalyst  layers  as  integral  parts 
of  the  unified  domains,  where  the  electrochemical  reactions 
are  incorporated  as  source  terms.  This  approach  eliminated 
the  necessity  of  prescribing  interfacial  boundary  conditions 
between  the  catalyst  layer  and  the  GDL.  Now  the  two-phase 
flow  model  is  also  applied  to  the  anode  side  and  the  two 
gas  channels.  In  addition,  the  momentum  transfer  between 
the  liquid  and  gas  phases  due  to  phase  change  is  taken  into 
consideration. 


2.  Experimental  system 

A  flow  diagram  of  the  experimental  system  is  shown  in 
Fig.  1.  The  fuel  cell  test  system  can  control  the  reactant  gas 
flow  rates,  fuel  cell  operating  temperature,  reactant  gas  hu¬ 
midification  temperature  and  operating  pressure.  The  reactant 
gas  flow  rates  are  regulated  by  two  mass  flow  controllers;  the 
humidification  temperatures  of  the  reactant  gases  are  con¬ 
trolled  by  regulating  the  humidifier  temperature;  and  the  op¬ 
erating  pressures  are  regulated  by  adjusting  the  back  pressure 
regulators.  A  nitrogen  purging  system  is  incorporated  in  the 
experimental  system  to  purge  both  the  anode  and  cathode 
sides  before  and  after  experiments  to  ensure  safety.  Experi¬ 
mental  data  are  recorded  after  the  fuel  cell  operation  reaches 
steady  state. 

In  this  study,  a  single  PEM  fuel  cell  with  an  active  area 
of  5  cm2  was  used.  The  membrane  used  was  Nation™  117 
membrane,  the  gas  diffusion  layers  (GDLs)  were  carbon  fiber 
paper,  and  the  carbon  supported  platinum  with  a  loading  of 
0.4  mg  cm-2  was  used  on  both  the  anode  and  cathode  sides. 
The  collector  plates  were  made  of  carbon  and  the  flow  fields 
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Computer  measure 
and  control  system 


Fig.  1.  The  schematic  of  the  fuel  cell  experimental  system. 


are  serpentine  on  both  sides.  The  end  plates  were  made  of 
copper.  Detailed  information  of  the  fuel  cell  used  is  listed  in 
Table  1. 


3.  Mathematical  model 

As  stated  earlier,  the  flow  model  used  is  the  mixture  two- 
phase  flow  model  [14,15],  and  we  used  the  two-phase  flow 
model  in  both  the  cathode  and  anode  sides,  including  the  two 
gas  channels.  The  catalyst  layers  are  also  included  in  the  uni¬ 
fied  domains  of  cathode  and  anode,  and  the  same  two-phase 
flow  model  is  applied.  In  addition,  the  momentum  transfer 
between  the  liquid  and  gas  phases  due  to  phase  change  is 
taken  into  consideration. 

The  modeling  domains  of  a  PEM  fuel  cell  are  shown  in 
Fig.  2.  Three  different  domains  are  considered,  one  for  the 
cathode  fluid,  one  for  the  anode  fluid,  and  one  for  water.  The 
cathode  and  anode  domains  each  consist  of  the  respective 
gas  channel,  GDL  and  catalyst  layer  (CL),  and  the  domain 


Table  1 

Geometric  parameters  and  catalyst  loading  of  the  experimental  fuel  cell 


Active  area  (cm2) 

5.0 

Channel  length  (cm) 

2.24 

Channel  width  (cm) 

0.1 

Channel  depth  (cm) 

0.09 

Gas  diffuser  thickness  (cm) 

0.03 

Catalyst  layer  thickness  (cm) 

0.00129 

Membrane  thickness  (cm) 

0.0108 

The  anode  and  cathode  CTL  loading  (mg) 

0.4 

for  water  include  all  the  components:  two  gas  channels,  two 
GDLs,  two  catalyst  layers  and  the  membrane. 

3.1.  Continuity  equation 

The  continuity  equation  for  the  two-phase  mixture  is: 

d(sp) 

+  V(spu)  =  0  (1) 

ot 

where  s  is  the  porosity  and  it  equals  one  in  the  channels  if 
the  liquid  volume  fraction  is  zero. 

3.2.  Momentum  equations 


The  momentum  equations  for  the  liquid  water  and  the  gas 
mixtures  in  both  anode  and  cathode  are  respectively: 


d(spu\) 

— - - b  S/{s\p\u\u\)  +  fn(u\ 

ot 


=  V(£l/ZlVwi)  + 


0  gas  channel 

-  —  (efm)  GDL  and  CL 
K\ 


d(spug) 


+  W  (SgPgUgUg)  +  fn(ug  —  u\) 


=  —SgS/P  +  V(£g/XgVwg) 


0 


gas  channel 


+ 


~  —  (slue)  GDL  and  CL 

K„  g  g 


where  fn  is  the  rate  of  phase  change. 
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Gas  channel  in  cathdoe 


Air  flow 
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Fig.  2.  The  unified  modeling  domains.  (1)  Domain  for  cathode  fluid;  (2)  domain  for  water;  (3)  domain  for  the  anode  fluid.  (L)  channel  length;  (D)  channel 
depth;  (< d)  gas  diffuser  thickness;  (5)  catalyst  layer  thickness;  ( h )  membrane  thickness. 


3.3.  Water  transport  equation  in  membrane 


Velocity  :  pu  =  p\u\  +  PgUg 


The  governing  equation  for  liquid  water  transfer  in  the 
membrane  is  given  by  [6] : 


Arnct  _  nd I  rA  dCw  n  Km  d Pw 

^vw  —  r-  i  cw  j 

F  dy  jl  dy 


3.4.  Species  conservation 


The  species  conservation  equation  is  given  by  [14,15]: 

d(spC)  rv 

+  V(yapwC  )  =  V(spDVCa) 
ot 


+  V 


sY/Pk^DtiyCl  -VCa)] 

k 


-  V 


Ec“^ 

k 

(5) 


where  k  denotes  either  gas  phase  or  liquid  phase,  a  indicates 
species,  such  as  oxygen,  nitrogen  or  water  in  the  cathode, 
and  it  indicates  hydrogen  or  water  in  the  anode,  and  ya  is  the 
correction  factor  to  be  given  below. 


3.5.  Mixture  properties  and  parameters 

The  properties  and  parameters  for  the  mixture  are  given 
below: 

Density  :  p  =  p\s  +  pg(l  —  s)  (6) 

Concentration  :  pCa  =  p\C^s  +  pgCg(l  —  s)  (7) 


Diffusion  coefficient  :  pDa  =  p\sD^  +  pg(l  —  s)D^ 


(9) 

p(hC«  +  XgC«) 


Correction  factor  :  ya  = 

'  D  6 

PisC[“  +  pg(l  -  s)C“ 

(10) 

Relative  mobility  :  X\  = 

krl/vi 

(ID 

kx\/v\  +  kxg/Vg 

Xg  =  \  —  X\  (12) 

1 


Fluid  kinematics  viscosity  :  v  =  -  (13) 

kri/vi  +  krg/vg 

Porosity  :  s  =  £\  +  £g  (14) 

where  s\  =  ss  representing  the  volume  fraction  of  liquid  wa¬ 
ter  in  the  porous  media,  Eg  =  (1  —  s)s  is  the  effective  poros¬ 
ity  that  denotes  the  volume  fraction  of  the  gas  mixture  in  the 
porous  media,  and  s  is  the  liquid  saturation,  representing  the 
volume  fraction  of  liquid  water  in  the  pores. 

In  addition  to  the  above  equations,  the  multiphase  mixture 
model  incorporates  some  important  relations  for  determining 
the  velocities  of  individual  phases  in  the  mixture  flow  field: 


sp\u\  =  N\  +  X\Epu 
spgUg  =  —N\  +  XgSpu 

X\Xp 

N\  =  K—^S/Pc 
v 


(15) 

(16) 

(17) 
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Ns  +  Ni  =  0 


3.6.  Constitutive  relations 


(18) 


The  porous  electrode  is  treated  as  homogeneous  porous 
media.  The  empirical  correlations  for  the  relative  permeabil¬ 
ities  for  the  liquid  and  gas  phases  are  respectively  [15]: 

=  ,3  (19) 

Kg  =  (1  -  s)3  (20) 

The  capillary  pressure  is  considered  as  a  function  of  satu¬ 
ration  and  the  empirical  equation  is  [14]: 

e  \  1/2 


where  R  is  oxygen  or  hydrogen;  k  equals  1  for  the  cathode 
side  and  0.5  for  the  anode  side. 

3. 7.  Relationships  for  the  membrane 

Water  may  transport  through  the  membrane  under  three 
different  transport  mechanisms:  electro-osmotic  drag,  back 
diffusion  and  hydraulic  permeation.  A  net  water  transport 
coefficient  is  defined  as  [e.g.  You  and  Liu] 


/  £  \ 

Pc  =  Pg-P\  =  <r(-J  m 

where  J(s\)  is  the  Leverett  function  given  by  [16]: 

J(s)  =  1.417(1  -  .v)  -  2.120(1  -  s)2  +  1.263(1  -  s)3 


(21) 


_  „  ,  A,  F  dCw  Km  F  dPw 

a  —  T  A  T  A 

22  I  ay  pi  I  ay 

The  water  diffusion  flux  in  membrane  is 
dC 


(27) 


A/'diff  —  _[) 


w 


(28) 


dym,w 

Water  molar  concentration  in  membrane  is  given  by  [18] 

Pm,dPw 7- 


Cw  = 


(22) 


(29) 


Pm.d  T  Pw  P\n 

Water  content  in  membrane  can  be  determined  by  water 
activity. 

xwP 


0.043  +  17.81a  —  39.85 a2  +  36.0a3  0 


<  a  = 


<  1 


sat 


A  =  14+  1.4 


22 


xwP 


sat 


(30) 


The  rate  of  phase  change  is  given  by  [17] 
D?ff  PMV  (  PK 


—  11  1  <  a  =  -  <  3 

/  ^sat 

X^P  Q 
a  =  -  >  3 

P sat 

The  water  diffusion  coefficient  in  the  membrane  can  be 
expressed  as  [18]: 


m  =  Sh 


—  /e - 

L  J  RT 


sat 


—  Vr 


(23) 


Dw  =  (11.102  +  9.129(1  —  exp  (— 0.108A.))) 
R  \  303  T 


(31) 


where  Sh  is  the  Sherwood  number,  xv  the  mole  fraction  of 
water  vapor,  P  the  pressure  of  the  gas  mixture,  and  Mv  the 
molecular  weight  of  water  vapor. 

The  saturation  vapor  pressure  is  given  by  [6] : 


-5^2 


log10  PSat  =  -2.1794  +  0.02953T  -  9.1837  x  10"JT 
+  1.4454  x  10_77’3 


The  ionic  conductivity  is  [18]: 

K  =  (-4.43  +  0.0179T)(— 0.0108  +  0.2365 
x  (1  -  exp  (-0.0285A))) 


(32) 


(24) 


The  liquid  water  saturation  s,  representing  the  volumetric 
fraction  of  void  space  occupied  by  the  liquid  water,  is  given 
by  [15]: 


The  membrane  resistance  is  determined  by  integrating  the 
resistance  over  the  membrane  thickness.  Taking  into  the  ef¬ 
fect  of  membrane  swell,  it  is  given  by: 


ARm  — 


s  = 


pC^°  -  pgcf° 

P\  -  PgCg2° 


J 


H(cdX  +  1)1/3  dyd 


o 


K 


(33) 


(25) 


where 


The  Butler- Volmer  equation  is  used  to  model  the  electro¬ 
chemical  reaction  rates  in  the  catalyst  layers 

k  _ 


Cd  = 


Pm.d^w 
Pw  Pm 


(34) 


3.8.  Boundary  conditions 


J  =  aio 


Cr 


C 


ref 

R 


exp 


RT 


exp 


—otc  Fr\ 
RT 


(26) 


Boundary  conditions  are  required  to  solve  the  continu¬ 
ity,  momentum  and  species  conservation  equations.  At  the 
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inlet  of  the  anode  and  the  cathode,  the  fluids  have  given  ve¬ 
locities.  The  standard  exit  boundary  and  no- slip  boundary 
conditions  are  used  at  the  channel  exits  and  channel  walls. 
In  the  species  field,  the  inlet  species  concentrations  are  pre¬ 
scribed,  the  species  gradients  are  set  to  zero  at  the  channel 
exits  and  at  the  walls.  The  gradients  of  reactant  (oxygen  or 
hydrogen)  concentrations  are  assumed  to  be  zero  at  the  inter¬ 
face  between  the  membrane  and  the  respective  catalyst  layer. 
The  detailed  boundary  conditions  are  given  below. 

At  the  cathode  inlet  (x  =  0) 


U  ljt=0  —  win> 


^lv=0 


=0  =  o,  cH2° 


x=0 


AbO  put 
^g.sat  R*bn 

Pg 


c° 2 


1O2 


x=0 


=  C~2,  0  <y<D 


where  RHj^  is  relative  humidify  of  inlet  air. 


^lr=0  0,  0 1  y=n  0, 


\x=0 
acH2° 

dx 


=  0, 


x=0 


dC°2 

dx 


x=0 


—  0,  D<y<D  +  d  +  8 

(36) 


where  RH“2  is  relative  humidify  of  the  anode  inlet  gas. 

acH2° 


^  I x=L  (t 

ac"2 


v\ x=L  =  0, 


dx 


=  0, 


x=L 


dx 


—  0,  D-\-d-\-8-\-h/£y/£D-\-  2  d  T  2  8  T  h 


x=L 


(42) 


a 

in 

du 

dx 

x=0 

(35) 

v\x= 

ii 

o 

II 

O 

At  the  anode  outlet  ( x  =  0) 

acH2° 


'dx 


=  0, 


X=0 


acHz 

dx 


X=0 


0, 

(43) 


u  |A=0 
dCu- 


acH2° 

dx 


=  0, 


x=0 


dx 


x=0 


—  0,  D  +  d  +  8  +  h<y<D  +  2  d  b  25  h 

(44) 


At  the  cathode  outlet  (x  =  L) 


du 

dx 


=  0, 


aCH2° 

dx 


x=L 

v\x=L  =°.  0  <  y  <  D 


=  0, 


x=L 


dc°2 

dx 


x=L 


=  o. 


(37) 


u\x=L=0,  V\X=L=0, 

dC°2 


dCn2o 

dx 


=  o, 


x=L 


dx 


—  0,  D^y^Df-df-S 


(38) 


x=L 


At  the  wall  of  the  cathode  channel  (y  =  0) 

^ly=0  ^|y=0  ^5 

acH2°  dc°2 

=  o,  - 

y=0  °y  y=0 


dy 


=  0,  0  <x<L  (39) 


At  the  interface  between  the  membrane  and  the  cathode 
catalyst  layer 

dc°2 


9 y 


=  0,  0  <x<L 


(40) 


y — D~\~d~\~S 


At  the  anode  inlet  (x  =  L) 

M  \j=L  ^H2,in’  v\x=l  0? 

] 

}g,sat 


ch2o 


pU2°.  rhh- 


m 


Cn2 


x=L  pg 

D  +  2d  +  28  +  h  <  y  <  H 


x=L 


=  c 


h2 


in 


(41) 


At  the  wall  of  the  anode  channel  (y-H ) 

3CH2° 


u\y—H  ~  0,  v\ y—H  —  0, 

acHz 


dy 


=  o, 


y=H 


dy 


=  0,  0  <x<L 


(45) 


y=H 


At  the  interface  between  the  membrane  and  the  anode 
catalyst  layer: 

acH2 


a  y 


=  0,  0  <x<L 


(46) 


y — D-\~d-\-8-\~h 


4.  Results  and  discussion 

The  governing  equations  with  the  boundary  conditions 
are  solved  numerically  by  a  finite-volume  method.  Table  2 

Table  2 

The  parameters  used  in  the  model  for  the  base  case 


Gas  diffuser  porosity 

0.4 

Gas  diffuser  tortuosity 

1.5 

Catalyst  layer  porosity 

0.25 

Air  stoic  ratio 

3.5 

Hydrogen  stoic  ratio 

2 

Catalyst  surface  area  (cm2  cm-1) 

1.4  x  105 

Exchange  current  density  (A  cm-2) 

4.84  x  10-8 

Cathodic  transfer  coefficient 

0.52 

Anodic  transfer  coefficient 

0.54 

Dry  membrane  density  (kgm-3) 

2.16 

Ionomer  equivalent  weight 

1100 

Sh  ft/L  (cm-2) 

20000 
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lists  the  parameters  used  in  the  base  case,  and  the  geometric 
parameters  are  listed  in  Table  1 . 

4.1.  The  effect  of  fuel  cell  operating  temperature 

It  is  well  know  that  operating  temperature  has  a  signifi¬ 
cant  effect  on  fuel  cell  performances.  The  activation  energy 
decreases  with  temperature  and  various  transport  properties, 
such  as  membrane  proton  resistance,  gas  diffusivities,  etc.  all 
depends  on  temperature.  Fig.  3  shows  the  modeling  results  of 
a  fuel  cell  performance  at  different  temperatures  compared 
with  experimental  data.  The  humidification  temperatures  on 
both  the  anode  and  cathode  sides  are  70  °C,  the  operating 
pressures  are  3  atm  on  both  sides,  and  the  flow  rates  are 
100  seem  on  the  anode  side  and  200  seem  on  the  cathode  side. 
Within  the  parameter  range,  the  performance  of  the  fuel  cell 
increases  with  temperature  and  the  modeling  results  agree 
well  with  experimental  data  except  at  high  current  density, 
where  the  model  over-predicts  fuel  cell  performance. 

Figs.  4  and  5  show  the  oxygen  mass  fraction  distribu¬ 
tions  in  the  cathode  catalyst  layer  at  the  same  overpotential 
for  two  different  operating  temperatures  of  70  and  80  °C,  re¬ 
spectively.  The  average  cell  current  densities  corresponding 
Figs.  4  and  5  are  0.926  and  1 .082  A  cm-2  respectively.  At  the 
same  overpotential,  the  cell  generates  higher  current  at  higher 
cell  temperature,  and  more  current  is  generated  in  the  region 
close  to  the  GDL  (low  value  of  y/5)  and  lower  oxygen  concen¬ 
tration  is  found  in  the  region  near  the  membrane  (high  value 
of  y/<5).  This  is  due  to  the  fact  that  at  higher  temperatures,  the 
chemical  reaction  occurs  more  readily  and  more  oxygen  is 
consumed  before  reaching  the  deeper  region  of  the  catalyst 
layer.  This  indicates  that  at  higher  operating  temperature  the 
effective  thickness  of  the  catalyst  layer  [13]  is  thinner  and 
this  result  should  be  taken  into  consideration  in  catalyst  layer 
design  and  optimization. 


Current  density  (A/cm2) 

Fig.  3.  Comparison  of  the  modeling  results  with  experimental  data  at 
different  fuel  cell  operating  temperatures.  Case  conditions:  H2  flow 
rate  =  100  seem,  humidification  temperature  on  both  sides  =  70  °C,  operat¬ 
ing  pressure  =  3  atm  on  both  sides;  air  flow  rate  =  200  seem. 


Fig.  4.  Oxygen  mass  fractions  in  the  cathode  catalyst  layer  at  fuel  cell  op¬ 
erating  temperature  of  70  °C.  Normalized  distance  is  from  the  cathode  cata¬ 
lyst  layer’s  interface  with  GDL  to  its  interface  with  the  membrane.  Case 
conditions:  H2  flow  rate  =  100  seem,  humidification  temperature  on  both 
sides  =  60  °C,  pressure  on  both  sides  =  3  atm;  air  flow  rate  =  200  seem. 

Figs.  6  and  7  show  the  total  water  mass  fraction  distribu¬ 
tions  across  the  cathode  catalyst  layer  at  different  locations 
along  the  flow  direction  for  two  different  operating  tempera¬ 
tures  of  70  and  80  °C,  respectively.  Total  water  mass  fraction 
increases  along  the  flow  direction  due  to  the  accumulation  of 
water  produced  in  the  catalyst  layer  and  the  electro-osmotic 
drag  from  anode.  Total  water  mass  fractions  are  generally 
higher  for  the  higher  temperature  case  (Fig.  7).  This  can  be 
attributed  to  the  higher  water  production  rate  and  higher  rate 
of  electro-osmosis  due  to  higher  current  densities. 


Fig.  5.  Oxygen  mass  fraction  in  the  cathode  catalyst  layer  for  fuel  cell  tem¬ 
perature  =  80  °C.  Normalized  distance  is  from  the  cathode  catalyst  layer’s 
interface  with  GDL  to  its  interface  with  the  membrane.  Case  conditions: 
H2  flow  rate  =100  seem,  humidification  temperature  on  both  sides  =  60  °C, 
pressure  on  both  sides  =  3  atm;  air  flow  rate  =  200  seem. 
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Fig.  6.  Total  water  mass  fractions  at  different  cross  sections  in  the  cathode 
catalyst  layer.  Normalized  distance  is  from  the  cathode  catalyst  layer’s  inter¬ 
face  with  GDL  to  its  interface  with  the  membrane.  Case  conditions:  H2  flow 
rate  =100  seem,  humidification  temperature  on  both  sides  =  60  °C,  pressure 
on  both  sides  =  3  atm;  air  flow  rate  =  200  seem.  Fuel  cell  temperature  =  70  °C. 

It  is  shown  in  Fig.  8  that  the  resistance  in  the  membrane 
for  the  fuel  cell  operating  temperature  at  80  °C  is  lower  than 
that  at  60  and  70  °C.  At  a  higher  operating  temperature,  the 
electrochemical  reaction  is  faster,  and  more  water  is  produced 
in  the  cathode.  More  water  makes  the  membrane  better  hy¬ 
drated,  and  thus  ionic  resistance  is  reduced.  Fig.  8  shows  that 
for  operating  temperature  at  60  °C,  resistance  in  the  mem¬ 
brane  decreases  along  the  flow  direction.  This  indicates  that 
there  is  more  water  in  membrane  near  the  outlet  than  that  near 
the  inlet. 


Fig.  7.  Total  water  mass  fractions  at  different  cross  sections  in  the  cathode 
catalyst  layer.  Normalized  distance  is  from  the  cathode  catalyst  layer’s  inter¬ 
face  with  GDL  to  its  interface  with  the  membrane.  Case  conditions:  H2  flow 
rate  =100  seem,  humidification  temperature  on  both  sides  =  60  °C,  pressure 
on  both  sides  =  3  atm;  air  flow  rate  =  200  seem,  fuel  cell  temperature  =  80  °C. 


Fig.  8.  Resistance  in  the  membrane  with  different  fuel  cell  operating  temper¬ 
atures.  X  is  the  flow  direction  and  X  =  0  is  the  inlet.  Case  conditions:  H2  flow 
rate  =100  seem,  humidification  temperature  on  both  sides  =  60  °C,  pressure 
on  both  sides  =  3  atm;  air  flow  rate  =  200  seem. 


Fig.  9  shows  water  content  distributions  across  the  mem¬ 
brane  at  different  locations  along  the  channel.  Water  content 
in  membrane  increases  from  anode  side  to  cathode  side  due 
to  electro-osmosis  and  water  generation  in  the  cathode.  It  is 
clear  from  Fig.  9  that  the  average  water  content  near  the  cath¬ 
ode  outlet  is  higher  than  that  near  the  inlet.  This  is  caused  by 
water  accumulation  along  the  channel  and  leads  to  a  lower 
ionic  resistance  near  the  outlet. 


Normalized  Distance 


Fig.  9.  Water  content  at  different  cross  sections  in  the  membrane.  Normal¬ 
ized  distance  is  from  the  membrane’s  interface  with  the  anode  CL  to  the  inter¬ 
face  with  the  cathode  CL.  Case  conditions:  H2  flow  rate  =100  seem,  humid¬ 
ification  temperature  on  both  sides  =  60  °C,  pressure  on  both  sides  =  3  atm; 
air  flow  rate  =  200  seem,  fuel  cell  temperature  =  60  °C. 
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Fig.  10.  Performance  comparisons  between  modeling  results  and  exper¬ 
imental  data  at  different  operating  pressures.  Case  conditions:  H2  flow 
rate  =100  seem,  humidification  temperature  on  both  sides  =  70  °C;  air  flow 
rate  =  200  seem,  fuel  cell  temperature  =  70  °C. 


4.2.  The  effect  of  operating  pressure 

Operating  pressure  is  also  an  important  parameter  for  PEM 
fuel  cell.  Fig.  10  shows  the  comparison  of  modeling  predic¬ 
tion  of  a  fuel  cell  performance  at  different  operating  pressures 
with  the  experimental  results.  The  fuel  cell  operating  temper¬ 
ature  and  humidification  temperature  on  both  the  anode  and 
cathode  sides  are  70  °C.  Inlet  flow  rates  are  100  seem  on  the 
anode  side  and  200  seem  on  the  cathode  side.  It  is  shown  in 
Fig.  10  that  modeling  results  agree  well  with  the  experimental 
results,  and  fuel  cell  performances  increases  with  operating 
pressure. 

Fig.  11  shows  the  oxygen  mass  fraction  in  the  cathode 
catalyst  layer  for  operating  pressure  at  2  atm.  A  comparison 
of  Figs.  5  and  1 1  clearly  shows  that  oxygen  mass  fraction  in 
Fig.  11  (lower  pressure)  is  higher  than  that  in  Fig.  5  (higher 
pressure),  and  the  gradient  of  oxygen  mass  fraction  is  lower 
in  the  lower  pressure  case  (Fig.  11).  At  a  higher  operating 
pressure,  the  electrochemical  reaction  rate  is  higher  and  the 
current  density  is  higher,  thus  more  oxygen  is  consumed  in 
cathode  catalyst  layer.  This  leads  to  lower  oxygen  concentra¬ 
tion  and  greater  oxygen  gradient  in  the  catalyst  layer. 

Fig.  12  shows  the  total  water  mass  fraction  at  different 
cross  sections  in  the  cathode  catalyst  layer  at  the  operating 
pressure  of  2  atm.  Comparing  Fig.  12  with  Fig.  7,  one  can 
see  that  total  water  mass  fraction  in  the  cathode  catalyst  layer 
for  a  lower  operating  pressure  (Fig.  12)  is  lower.  This  can  be 
attributed  to  the  lower  water  generation  rate  and  the  lower 
electro-osmosis  due  to  the  lower  current  density. 

Fig.  13  shows  the  variation  of  liquid  water  saturation 
across  the  cathode  channel  and  GDF  at  different  operating 
pressures.  It  is  clear  that  liquid  water  saturation  increases 
with  operating  pressure.  Fiquid  water  appears  when  the  wa¬ 
ter  vapor  partial  pressure  exceeds  the  local  water  saturation 


Fig.  11.  Oxygen  mass  fractions  in  the  cathode  catalyst  layer  with  operating 
pressure  at  2  atm.  Normalized  distance  is  from  the  cathode  catalyst  layer’s 
interface  with  GDL  to  the  interface  with  the  membrane.  Case  conditions: 
H2  flow  rate  =100  seem,  humidification  temperature  on  both  sides  =  60  °C, 
pressure  on  both  sides  =  2  atm;  air  flow  rate  =  200  seem,  fuel  cell  tempera¬ 
ture  =  80°  C. 

pressure.  At  a  higher  pressure,  the  vapor  faction  in  the  gas 
stream  is  lower  and  more  water  exists  in  liquid  form.  Besides, 
at  a  higher  operating  pressure,  the  fuel  cell  current  density  is 
higher  and  more  water  is  produced.  It  is  also  shows  in  Fig.  13 
that  the  liquid  water  saturation  in  the  GDF  is  larger  than  that 
in  the  gas  channel. 

4.3.  The  effect  of  humidification  temperature 

Humidity  of  input  gas  is  another  very  important  pa¬ 
rameter  for  PEM  fuel  cell  operations.  For  most  humid- 


Normalized  Distance 


Fig.  12.  Total  water  mass  fractions  at  different  cross  sections  in  the  cathode 
catalyst  layer.  Normalized  distance  is  from  the  cathode  catalyst  layer’s  inter¬ 
face  with  GDL  to  the  interface  with  the  membrane.  Case  conditions:  H2  flow 
rate  =100  seem,  humidification  temperature  on  both  sides  =  60  °C,  pressure 
on  both  sides  =  2  atm;  air  flow  rate  =  200  seem,  fuel  cell  temperature  =  80  °C. 
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Fig.  13.  Liquid  saturation  at  x  =  22.36  mm  in  the  cathode  with  different  op¬ 
erating  pressures.  The  coordinate  y  starts  at  the  channel  wall  (see  Fig.  2). 
Case  conditions:  FG  flow  rate  =  100  seem,  humidification  temperature  on 
both  sides  =  80  °C;  air  flow  rate  =  200  seem,  fuel  cell  temperature  =  80  °C. 

ifiers,  humidity  is  regulated  by  controlling  the  temperat¬ 
ure. 

Fig.  14  shows  the  liquid  water  saturation  in  the  cathode 
at  different  humidification  temperatures.  The  humidification 
temperatures  on  the  anode  side  for  each  case  equal  that  on 
the  cathode  side.  It  is  seen  from  Fig.  14  that  the  liquid  water 
saturation  increases  with  the  gas  stream  humidification  tem¬ 
perature.  This  is  an  obvious  result,  since  more  water  vapor 
enters  the  fuel  cell  at  higher  humidification  temperatures. 

Fig.  15  shows  the  total  water  mass  fraction  at  the  interface 
between  the  anode  GDL  catalyst  layer  at  three  different 
humidification  temperatures.  Note  that  the  anode  flow  is 
in  the  opposite  direction  to  that  of  the  cathode  due  to  the 
counter-flow  arrangement.  Under  this  low  hydrogen  flow 
rate,  along  the  anode  flow  direction,  total  water  mass  fraction 


Y  (mm) 

Fig.  14.  Liquid  saturation  atx  =  22.36  mm  in  the  cathode  channel  and  GDL  at 
different  humidification  temperatures.  The  coordinate  y  starts  at  the  channel 
wall  (see  Fig.  2).  Case  conditions:  H2  flow  rate  =  100  seem,  pressure  on  both 
sides  =  3;  air  flow  rate  =  200  seem,  fuel  cell  temperature  =  70  °C. 


X  (mm) 

Fig.  15.  Total  water  mass  fractions  at  the  interface  of  GDL  and  CL  on  the  an¬ 
ode  at  different  humidification  temperatures.  X  represents  the  flow  direction 
and  x  =  0  is  the  inlet  (see  Fig.  2).  Case  conditions:  H2  flow  rate  =  30  seem, 
pressure  =  3  atm  on  both  sides;  air  flow  rate=  120  seem,  fuel  cell  tempera¬ 
ture  =80  °C. 


X  (mm) 

Fig.  16.  Liquid  saturation  at  the  interface  of  GDL  and  CL  at  the  anode  at 
different  humidification  temperatures.  The  coordinate  y  starts  at  the  channel 
wall  (see  Fig.  2).  Case  conditions:  H2  flow  rate  =  30  seem,  pressure  =  3  atm 
on  the  both  sides;  air  flow  rate  =120  seem,  fuel  cell  temperature  =  80  °C. 

increases.  This  is  due  to  the  fact  that,  for  each  mole  of 
hydrogen  consumed,  the  net  water  transfer  to  the  cathode 
side  is  less  than  one  mole.  It  can  also  be  seen  from  Fig.  15 
that,  the  total  water  mass  fraction  is  higher  at  higher  gas 
humidification  temperatures. 

Fig.  16  shows  the  liquid  water  saturation  in  the  anode  cat¬ 
alyst  layer  at  its  interface  with  GDL.  Liquid  water  increases 
along  the  anode  flow  direction  due  to  hydrogen  consumption 
as  discussed  above.  Again,  at  higher  humidification  temper¬ 
atures,  liquid  water  saturation  is  higher. 

5.  Conclusion 

A  two-phase  flow  PEM  fuel  cell  model  is  developed  for 
the  PEM  fuel  cell.  The  mixture  two-phase  flow  model  is  used 
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in  both  the  cathode  and  the  anode,  and  it  also  used  in  the  two 
gas  channels  in  addition  to  the  two  GDLs  and  two  catalyst 
layers.  The  modeling  results  at  different  operating  conditions 
are  compared  with  experimental  data  and  good  agreements 
are  obtained.  The  model  is  used  to  study  the  effects  of  vari¬ 
ous  operation  parameters,  and  specifically  the  results  of  water 
distribution  across  different  parts  of  a  fuel  cell  are  presented. 
The  following  conclusions  can  be  obtained  from  this  study, 
(a)  In  the  cathode  catalyst  layer,  total  water  mass  fraction  in¬ 
creases  from  its  interface  with  the  DGL  to  its  interface  with 
the  catalyst  layer.  Within  certain  temperature  range  and  suf¬ 
ficient  gas  stream  humidification,  total  water  mass  fraction 
in  the  cathode  catalyst  layer  and  the  membrane  ionic  con¬ 
ductance  increases  with  fuel  cell  operating  temperature  due 
to  the  higher  water  generation  rate,  (b)  The  total  water  mass 
fraction  in  the  cathode  catalyst  layer  and  the  liquid  volumet¬ 
ric  fraction  (liquid  saturation)  in  the  cathode  channel  and 
DGL  increase  with  operating  pressure  due  to  higher  water 
generation  rate  and  lower  vapor  faction  in  the  gas  phase,  (c) 
For  counter-flow  arrangement,  generally,  total  mass  fraction 
and  liquid  water  volumetric  fraction  in  the  anode  GDL  and 
catalyst  layer  increases  along  the  anode  flow  direction. 
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